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The role of the electron-hole symmetry breaking is investigated for a symmetrical commutative 
two-level system in a metal using the multiplicative renormalization group in a straightforward way. 
The role of the symmetries of the model and the path integral technique are also discussed in detail. 
It is shown that the electron-hole symmetry breaking may make the model non-commutative and 
generate the assisted tunneling process which is, however, too small itself to drive the system into 
the vicinity of the two-channel Kondo fixed point. While these results are in qualitative agreement 
with those of Moustakas and Fisher (Phys. Rev. B 51, 6908 (1995), ibid 53, 4300 (1996)) the scaling 
equations turn out to be essentially different. We show that the main reason for this difference is 
that the procedure for the elimination of the high energy degrees of freedom used by Moustakas and 
Fisher leaves only the free energy invariant, however, the couplings generated are not connected to 
the dynamical properties in a straightforward way and should be interpreted with care. These latter 
results might have important consequences in other cases where the path integral technique is used 
to produce the scaling equations and calculate physical quantities. 
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I. INTRODUCTION 



The two-channel Kondo model originally introduced for magnetic ions with crystalline splitting)!! and tunneling 
two-level systems (TLS) in metalaJ has recently attracted broad interest.Q'El In some of those problems the orbital 
variables play the role of the spins of the impurity and the conduction electron while the real spins of the conduction 
electrons are not dynamical variables but instead represent an extra twofold degeneracy for the conduction electrons 
interactiiig|-.with the impurity. A TLS can be considered as a good realization of the two-channel orbital Kondo 
problemBt! Besides the fast TLS's and the impurity ions with orbital dynamicaLyariables mentioned before the two- 
channel Kondo problem can also be related to the two-impurity Kondo problem. alilll Due to the double degeneracy of 
the real conduction electron spins the fixed point of these models is in the intermediate strong coupling region.^ 

It became obvious in recent years that the breaking-of-electron-hole (e-h) symmetry may result in the instability 
of the fixed points found in ±he e-h symmetric case,II3II3 and new fixed points might become relevant as in the 
two-impurity Kondo problemflo 

The -present paper is devoted to a special class of the TLS problem, called the dissipative or commutative TLS 
modelJiilO where there is only one dynamical coupling to describe the electron-TLS interaction in contrary to the 
orbital Kondo problem (non-commutative TLS problem), where there must be at least two non-commuting interaction 
terms. Usually, when dealing with such impurity models an e-h symmetrical conduction electron band is assumed. 
Then the dissipative two-state model provides an example of marginal models in the renormalization group sense, 
i.e., its dynamical couplings remain unrenormalized under scaling. The main goal of the present paper is to show 
that the commutative model becomes unstable due to the presence of e-h symmetry breaking of a specific form, and 
serves as a toy model to examine the role of e-h symmetry breaking occurring in a realistic model of a TLS. As we 
shall see, in the presence of e-h symmetry the model has an additional symmetry which ensures the stability of the 
line of marginal fixed points. This problem has its own interest even if it does not play a major physical role in most 
of the cases of realistic systems. By developing numerical estimates, we will argue that for realistic parameters of a 
TLS, this effect is extremely small for metals, and can be safely neglected for practical purposes. Therefore, in case 
the screening interaction is not anomalously strong, the dominant mechanism making the commutative TLS model 
unstable will be the intrinsic assisted tunneling considered in Refs. [| and 

The essential effect of the e-h symmetry breaking can be demonstrated by considering the two time-ordered diagrams 
in Fig. [I]. If the conduction electron is scattered by a simple potential scatterer V, both diagrams provide a diverging 
term ~ V 2 ln-^j, (D and u> being the conduction electron bandwidth and some small-energy scale, respectively), but 
these contributions cancel each other due to e-h symmetry.O That is one of the most important ingredients of the 
X-ray absorbtion problem. t3 In case of e-h symmetry breaking the cancellation is not exact, and a non-diverging 
contribution remains. Such contribution may change the mathematical structure and the universality class of the 
theory. 

To be specific, we imagine a TLS, where a heavy particle (HP) can jump between two positions in a double potential 
well. In general the HP can dynamically interact with the conduction electrons in two different ways (see Fig. ||): 

(i) The electrons try to form a screening cloud around the heavy particle. If only this interaction is taken into account 
then, using the logarithmic approximation, the usual e-h symmetrical modela (called the dissipative or commutative 
TLS model mentioned above) results in a marginal theory, where the screening interaction is unrenormalized but the 
hopping rate of the HP is essentially reduced, because the overlap matrix element of the electronic screening clouds 
at neighboring positions of the HP vanishes £3 

(ii) The conduction electrons also assist the-hopping of the HP (assisted hopping) both in the case of tunneling and 
virtual transitions through the excited states £3 If the second process (ii) is also taken into account then the couplings 
are not commuting in the momentum space, and the non-commutative model belonging to the class of two-channel 
Kondo models is recovered. 

In the present paper we concentrate on case (i), where the assisted tunneling (ii) is absent in the bare hamiltonian. 
As mentioned above, in the original treatment of the problem one usually assumes that the local density of states 
is the same for each orbital channel of the conduction electrons and that it is e-h symmetrical. In a realistic case, 
however, these assumptions are not a priori justified. i— in 

Recently, using a path integral approach Moustakas and FisherEirE3 argued that in the commutative problem (i) 
the above mentioned potential scattering at the HP is relevant and changes the scaling behavior of the TLS. The 
possibility has been pointed out that such a potential may generate an electron assisted hopping. Such a potential is 
usually neglected which is certainly not justified for the case when the heavy particle is different from the atoms of the 
host material. (The problem of this static potential scattering has been investigated earlier by Kagan and Prokof 'ev 
by applying the adiabatic renormalization technique.t£l) 

As we shall see both scattering on a local static potential and a realistic dispersion relation for the conduction 
electrons validate the previous assumptions for the local density of states, and generate process (ii). Therefore, in the 
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present paper we shall drop them and treat a general local density of states incorporatring both mechanisms. 

To investigate the relevance of the static potential and the generation of assisted tunneling two different methods 
will be developed. First, applying a straightforward multiplicative renormalization group method in the leading and 
next to leading logarithmic approximation, appropriate in the small coupling limit, we show that while the energy 
dependence of the local electronic spectral functions really generates the assisted tunneling term, it can be neglected 
in most of the cases and the commutative behavior is recovered. We also investigate the problem with a path integral 
technique. In the concluding section of this paper and Appendix B we discuss the ambiguity in the construction 
of the scaling equations derived by Moustakas and Fisher. We point out that in the path integral approach the 
elimination of the events with small time difference is not uniquely defined, and two procedures are possible: 6Q Any 
events being closer to each other than a certain time difference must be replaced by a generated interactions, (ii) 
Only those processes must be eliminated, where electrons are involved. While both procedures leave the free energy 
of the HP unchanged, the couplings generated and the electronic sector behave essentially differently. We will argue 
that the physical argument behind the second procedure is more sound. 

It is interesting to note that the first method applied by Moustakas and Fisher generates the electron-TLS coupling 
even in that trivial case where there is no dynamical coupling between the TLS and the conduction electrons.EUEa (This 
importaint point has also been mentioned in footnote 20 of Ref. p2| ) . We show that the meaning and the application 
of these couplings to calculate real physical quantities is not straightforward, as the scaled couplings are usually 
plugged into formulas obtained by conventional diagram technique not directly related to the path integral method. 
In this way different sets of the scaled couplings lead to different measurable quantities. Thus in the path integral 
formalism no more than one if any scaling procedure can be appropriate to provide the correct scaled couplings to 
calculate physical quantities. Since the path integral technique should be also correct in the weak coupling region 
where the multiplicative renormalization group technique is adequate, we can decide which procedure is correct by 
comparing the scaling equations in this limit. We show that the second scaling procedure is the one which reproduces 
the dynamical behavior obtained in the next to leading logarithmic approximation. With the elimination procedure 
(ii) the assisted hopping is not generated by the combination of spontaneous tunneling and simple potential scattering 
if all the dynamical couplings are switched off. 

We also show using the multiplicative renormalization group method that in the presence of e-h symmetry and 
without assisted tunneling the TLS-conduction electron dynamical couplings remain unrcnormalizcd up to the next 
to leading logarithmic order. This also contradicts to the elinimation scheme used in Refs. pT| , p2"[ 

The relevance of .different terms in the Hamiltonian has been discussed by Moustakas and Fisher in terms of 
symmetry breaking rWA l n order to clarify the problem special emphasis will be put on these symmetries. We show 
that the conventional commutative TLS model without potential scattering and with an e-h symmetrical band has an 
additional symmetry reflecting the role of e-h symmetry. This symmetry is a combination of reflection to the middle 
of the TLS for the conduction electrons (but not for the tunneling HP), and the e-h symmetry, and it can be broken 
in different ways, e.g. by potential scattering, an asymmetrical conduction electron band, or an energy dependent 
coupling between the HP and the conduction electrons. (This latter, analized in Appendix B of Ref. |l| is closely 
related to the energy dependence of the local density of states, studied in the present paper.) in the Hamiltonian which 
do not conserve the number of right (left) electrons in the sense introduced by Moustakas and Fisher O The breaking 
of e-h symmetry can a ls o, b e the consequence of an energy dependent coupling in agreement with the suggestion of 
Moustakas and FisherElH The new symmetry found guarantees that the low energy fixed point is marginal (the 
Hamiltonian is invariant under the scaling procedure) and no Kondo effect occurs. On the other hand if it is violated 
then in general the assisted tunneling is generated and the system starts to flow in the direction of the two-channel 
Kondo fixed point until the energy scale becomes smaller than the renormalized splitting between the two lowest 
states of the TLS, and the TLS dynamics is frozen out. 

However, we shall see that the assisted tunneling generated in this way is very weak and usually unobservable. In 
fact, the numerical renormalization group study of the same problem by Libero andpOliveira in the absence of e-h 
symmetry does not indicate the scaling toward the intermediate coupling fixed point Ba This is due to the fact that 
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since the generated assisted tunneling is always very weak the splitting of the TLS stops the scaling before the system 
could approach the neighborhood of the two-channel Kondo fixed point. It will be speculated, that it can result in an 
observable effect only in the case of very asymmetrical band (degenerate semiconductors )£-3 

In general the question can be raised whether a spontaneous hopping process renormalized by the screening interac- 
tion between the TLS and the electrons can induce an effective assisted hopping or not. A straightforward calculation 
of the diagrams in Fig. || shows that electron assisted hopping is generated but only if the local electronic propagator 
breaks e-h symmetry The simplest realization of such symmetry breaking is an asymmetric electron band or a static 
potential scattering at the tunneling centers. The diagram in Fig. [| is logarithmic, but it contains an extra small 
factor ~ Ao/ef ~ 10 -4 , A being the amplitude of spontaneous tunneling of the TLS and €f the Fermi energy. Such 
processes are usually dropped in the next to leading logarithmic calculations because of the appearance of the small 
term A /D in these equations. The appearance of this instability can be traced back in the path integral technique 
and supports the legitimity of the the second elimination scheme discussed above. 

The paper is organized as follows. In Sec. II the model is introduced and the time-ordered perturbation theory is 
applied to demonstrate the importance of the e-h symmetry breaking and how the assisted tunneling is generated. 
In Sec. Ill the additional symmetry is introduced in the presence of e-h symmetry and it is shown how that insures 
the stability of the marginal fixed point. In Sec. IV the strength of the assisted tunneling generated by the breaking 
of the e-h and the additional symmetry is estimated. In Sec. V the results and the ambiguity in the path integral 
method are discussed. In Appendices A and B we calculate the amplitude of the different e-h symmetry breaking 
parameters and e show how the scaling equations obtained in the different path integral techniques may be related 
by a formal transformation. Finally, in Appendix C and D we investigate the mixing of the high- and low-energy 
conduction electron degrees of freedom due to a simple potential scatterer and the e-h symmetry in the case of a cubic 
tight binding model. 



II. THE MODEL AND APPLICATION OF THE TIME ORDERED DIAGRAM TECHNIQUE 



First we introduce the model and the notations used closely following Ref. 22 . Since the results of this Section are 
independent of the real spin indices of the electrons, for the sake of simplicity we drop them throughout this Section. 
The creation operators of the HP at site 1 and 2 are denoted by df and ■ The conduction electron annihilation 
operators at positions ±R/2, c±, can be expressed in terms of the annihilation operators for spherical waves (s- waves) 
with momentum k at the positions r = ±R/2 denoted by c+,fc and C-^: 

f dk f dk f k 2 di}u , . „ , . , 

c ± = j -c ±k = J^J^ ¥ exp( ±i kR/2 )ck , (2.1) 

where the Ck's denote plane wave annihilation operators. Furthermore for convenience the following even and odd 
operators are defined: 



c 



• '■• — / at n \ ( C + 



(<*+<*)> 



where the normalization constants are given by 



c+ =__( C + -c+ ), (2.2) 



, r 2k 2 / sinkR\ 

N eA k) = - 1 + -^] , (2-3) 
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and the operators in Eq. (2.2) satisfy the anticommutation relations: {c ak , c^,} — 2n5(k — k')5 a p with a, (3 = e,o. 
The spin variables are not indicated in the above formulae. Then the Hamiltonian can be written in a concise form 
in terms of the operators 



e(o) 



dk 



(2.4) 



2?r e(o ^ k ' 

Assuming a local HP-electron interaction, the Hamiltonian considered by Moustakas and Fisher can be written as£H£3 

H = H Q + H d + U , (2.5) 



where 



dk 
2^ 



e fc (c+ fc c efc + c+ k c ok ) 



H d = A (dtd 2 +dtd 1 ) , 
U = Vi(c+c e + c+c ) + V 2 (ctc e - c+c D ) + V 3 (d+ di ~ d+d 2 ){ctc + c+c e ) 



(2.6) 

(2.7) 
(2.8) 



Here e k is the energy of the conduction electrons, Ao is the spontaneous tunneling rate of the HP, and the different 
momentum dependent couplings have been approximated by their values at the Fermi energy. The terms proportional 
to V\ and V% describe potential scatterings while V3 denotes the amplitude of the screening interaction. Note that 
in the interaction part U the constraint dfdi + d\d 2 = 1 has already been taken into account. In the case of an 
asymmetric double potential well there is an additional asymmetry term for the HP 



H' d = A z (d+d 1 -d+d 2 ) . 



(2.9) 



The couplings V\ and V 2 can be incorporated to the one-particle properties of the conduction electrons. The 
unperturbed one particle Green's functions are usually defined as 



Gi° fe ) K)=G^K) = 



(o). 



1 



(2.10) 



where the ui n 's denote Matsubara frequencies. Then the electronic local Green's functions modified by the potential 
scattering 



Sap G a (uj n ) = -{T{c a c~p}) u 
with a — e, o can be calculated easily and one obtains 



1 - (Vi ± V 2 )G$l(u n ) 



(2.11) 



(2.12) 



where ci°o(^) — J dk (iuj — e(fc)) 1 /27r. The corresponding spectral functions are given by 



f?e,o 



TmG$l(uj-iS) 



l-(l/i±l/ 2 )ReGi ](^-z<5) 



(Vi±V 2 )ImG^(w-t*) 



(2.13) 
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Since to the lowest order in uj the imaginary and real parts of the Green's functions can be approximated as ImGt?l(u) — 
iS) ~ const and Red ] (a; — iS) ~ u>, from Eq. ( 2.13 ) immediately follows that in the presence of V% or V% the e-h 



symmetry can not hold even if the original band is symmetrical. 
Now it is useful to introduce the Pauli operators 

d+d 2 - d+d! = t z , (2.14) 
dtd 1 +d+d 2 =T x (2.15) 

to describe the motion of the TLS, then the remaining terms to be treated are 

V 3 (e+c + c+c e )T z + A t x . (2.16) 

To make the analogy with the two-channel Kondo model more transparent a new rotated representation of spin 
operators can be introduced, t x — > f z and r z — > —f x and then the interaction Hamiltonian becomes 

- V 3 (ctc + c+c e )f x + A f z . (2.17) 

That interaction has resemblance to the anisotropical Kondo Hamiltonian in an external field Ao with a single coupling 
J x , but with Q e (u) ^ Q (u). 

As we have shown above the effect of the the scattering processes V\ and V2 can be incorporated into the spectral 
functions g e ^ {uj), so in the following we formulate everything in terms of these quantities and drop the terms V\ 
and V2. Since these are slowly varying functions, as a first step, they can be approximated by their values at the 
Fermi energy. Then treating V3 perturbatively one can notice that in the different diagrams only the combination 
(Pe(0)£> o (0)) V 2 occurs and therefore the corresponding dimensionless coupling is given by 

9=(Qe(0)g o (0)) 1/2 V 3 . (2.18) 

Apart from the present modification in this approximation the calculation is identical to the one in Ref. |l4| Calculating 
the leading and the next to leading logarithmic diagrams no corrections to the invariant couplings occur as the second 
order vertex corrections cancel the self-energy correction part. The only relevant renormalization is due to the HP wave 
function renormalization which generates the renormalization of the spontaneous splitting: Aq(uj) — Aq(cj/D) Nb9 / 2 , 
where D is a high energy cutoff of the order of the Fermi energy, ep, and N s = 2 denotes the spin degeneracy of the 
conduction electrons O 

The next question is whether the uj dependence of g e ^ (uj) due to the breaking of e-h symmetry is crucial in these 
considerations or not. As it can be seen from Fig. [l]the diagrams for first order vertex corrections are diagonal both 
in the even-odd electron channels and the TLS indices, and therefore they only renormalize and f> e ( )(w). Their 
contribution can be estimated by approximating 

Pe,oH = 0e,o(O)(l +a e ^uj) , (2.19) 

where a e ^ a a « 0.1 e^ 1 — 0.3 e^ 1 . Then the logarithmic part of the diagrams in Fig. |l| is exactly canceled and they 
give only a constant renormalization of the potential scattering, which can be taken into account by the redefinition 
of the a ei0 's.c3 The order of magnitude of this renormalization can be estimated as 

8{g{Q)V ia ) - g 2 aD < 1 . (2.20) 

It is important to note that the constants a e ,o characterize the spectrum of the local low energy excitations coupled 
to the TLS, and they should remain unsealed under the elimination of the high energy states D — > D' . This subtility 
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turns out to be non-trivial if one tries to simulate the e-h symmetry breaking with a potential scattering term as in 
Refs. |2l|j22| since then a ~ 1/D (see Appendix C). The main difficulty arises from the fact that the potential scattering 
strongly mixes the original low- and high-energy excitations. One could, in principle, proceed in two different ways: (a) 
One tries to treat the potential scattering similarly to the dynamical interaction terms and eliminates the high-energy 
degrees of the model without potential scattering (D — > D'). This procedure has been used in Refs. |2l],|2^. Since 
a ~ 1/D this procedure obviously changes the local low-energy excitation spectrum coupled to the TLS. (b) One first 
diagonalizes the exactly solvable potential scattering part of the Hamiltonian to find the real low- and high-energy 
excitations coupled to the TLS and then eliminates the latter ones. This is the one followed throughout this paper. 
Some insight to the difference between these two procedures is given in Appendix C. As shown in Appendix C, due 
to the strong mixing of the original low-energy and high-energy degrees of freedom, procedure (a) also eliminates a 
fraction of the real low energy excitations (in the presence of potential scattering) coupled to the TLS, and therefore 
the correct procedure is the second one. The essential difference between the two elimination procedures is closely 
related to the fact that the mass terms (i.e. the potential scattering in our case) can never be treated at equal footing 
with the interaction terms in a renormalizable field theory.Ej We also mention that elimination scheme (a) results in 
scaling equations which do not reproduce the pcrturbativc results. 

Let us now observe that the dia grams in Fig. 3 generate a new assisted tunneling process (ii) originally absent from 
the starting Hamiltonian Eq. (21): 



iJ asst . = iAir x (c+c e - c+c ) = -Air a: (cfc2 + C+Ci) , 



(2.21) 



where the annihilation operators for the orthogonal left and right electron states ci 2 have been defined following 
Ref. E2I as 



Cl,2 = (c e ± C )/V2 



(2.22) 



The contributions of the diagrams in Fig. 3 to the scaling equations depend on whether the intermediate electron line 
is even or odd and can be calculated by rescaling the cutoff D — ► D — dD (see Ref. ||) 



If - 


p-D+dD' 


J D — dD 


l-D 



1 



£e !O (0)(l + a e , e)de , 



(2.23) 



and generate the following terms 

.dD 



D 



g e a e 



Qo^o { 4. — I — \ Qe^-e 
( c e c e + CjC ) H 



{4 



(2.24) 



Here the first term is a dynamical renormalization of the tunneling amplitude which commutes with the screening 
term V3, but the second term is a generated assisted tunneling. 

Thus if g e a e ^ g a then the two diagrams in Fig. 3 generate the assisted tunneling process. For g e — g Q the 
condition a e 7^ a is equivalent with the one that in the site representation cx,2 the Green's function for the conduction 
electrons has also offdiagonal matrix elements. This result is in accordance with the suggestion of Moustakas and 
Fisher who related the appearance of the assisted tunneling to cr oss sc attering between the electron channels T' an d '2' 
and in this way to the non zero value of the coupling V2 (see Eq. ( 2.13| )).Ell However, as one can see from Eq. (2.24) one 
has a logarithmic contribution only if the electron-hole symmetry is broken in the even and odd channels in different 
ways, which means that this is a combination of 1 — > 2 scattering with the e-h symmetry breaking which generates 
the assisted tunneling. (In the simple model of Ref. |2l| the term V2 generated both cross scattering between channels 
1 and 2 and the e-h symmetry breaking.) Another important difference is that the term generated is proportional to 
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g 2 in contrast to the scaling equations of Ref. ^l],^2|, where such a coupling is also generated in the 5 = case, where 
the conduction electron sector and the HP is dynamically completely decoupled. This indicates that, although the 
scaling equations of Ref. ^l] leave the free energy of the system invariant, the scaled couplings there have nothing to 
do with physically measurable quantities like the scattering amplitudes directly connected to the vertex function, and 
must be interpreted with care. _ 
Thus for a e ^ a the assisted tunneling term is really generated and it makes the marginal fixed point unstable.EI 
However, its coefficient is very small as g ~ 0.1 and Aq/Dq ~ 10~ 4 , and therefore it has always been dropped in 
the previous calculations, even if it gives a logarithmic contribution (see Refs. |l4],^j). It will be shown later that this 
amplitude is too small to generate an assisted tunneling which brings the TLS into the vicinity of the two-channel 
Kondo fixed point described in the introduction when the hopping rate Ao is also taken into account as a lower cutoff 
for the scaling. 

It is worth mentioning that the c-h symmetry can also be broken in the case of a constant density of states but by 
using different upper and lower cutoffs D upp and — D\ ow . In this case, however, the contribution is not logarithmic. 



III. ANALYSIS OF THE ADDITIONAL SYMMETRY VALID FOR ELECTRON-HOLE SYMMETRY IN 
THE ABSENCE OF THE ELECTRON ASSISTED HOPPING Ai 

In the present section we analyze the symmetry properties of the commutative TLS problem and show that the 
model has an additional very restrictive symmetry related to the generation of the assisted tunneling. This symmetry 
is closely related to the-4j»*iscrvation of the number of the conduction electrons in channels '1' and '2' introduced by 
Moustakas and Fisher .BB 

The formulation of the e-h symmetry needs special care. In the general case the e-h symmetry does not hold 
exactly. There is, however, an approximate transformation nearby the Fermi surface which connects electron and hole 
states at the same part of the Fermi surface. In this approximate transformation the density of states is assumed 
to be independent of the energy and the momentum dependence of the different phase factors and wave functions is 
neglected, therefore it holds only in a very restrictive sence. The transformation is illustrated in Fig. ^, where a small 
section of the Fermi surface at the Fermi wave vector is considered. 

The deviations in the normal directions at that small section are neglected and the dispersion is linearized along 
the normal direction. Thus the electron state with momentum k and energy -E(k) can be characterized by a wave 
vector at the Fermi surface, kp and the energy e = E(k) — Ep: 

k = k F + VF(kF) e, (3.1) 

where v^(k^) is the Fermi velocity at point hp. Then the e-h transformation can be written as 

c kF , e -> 4 F _ £ . (3.2) 
In that approximation the energy independent density of states can be written as 

go = ji^J\^pj\> (3 - 3) 

where dSk F is the Fermi surface element. That approximation becomes more accurate as the Fermi surface is ap- 
proached. Thus a systematic expansion around the Fermi surface can be developed, where the zero order contribution 
is given above and the further terms can be given in terms of e. Such expansion has been used in the previous Section 
for the local density of states by introducing the a e o factors. 



As an illustration the electron annihilation operators at ±R/2 are c onsi dered for |Rk,p| <C 1, thus the phase factors 
are slowly varying functions of e. Now c± is defined similarly to Eq. ( |2.l|) : 



c± = I P- f ^|^ e ±l(k ^ ev(kF)/v2(k ^ )R/2 c kF , e . (3.4) 
J 2n J {2ny v[k F ) 

Ignoring th e e- depen dencc of the exponent the e-h stransformation takes the simple form c± — > (czp) + , which combined 
with Eqs. (2.2) and (2.22) gives the simple form 

ci — > . (3-5) 

In the calculation of the infrared divergencies in the leading order the approximate e-h symmetry can be applied. The 
corrections in the sense discussed above are less divergent or even convergent due to the occurring extra powers of 
e. Thus the above approach provides a systematic expansion in breaking of the e-h symmetry for an arbitrary Fermi 
surface. 

There are, however, special cpsf*i where exact e-h transformations exist, e.g., for a half filled cubic tight binding 
model discussed in Appendix 0113 til These transformations differ essentially from the above-discussed one as different 
regions of the Brillouin zone are connected by them. 

In order to study the interaction first its structure must be discussed. The following combinations of the creation 
and annihilation operators of the HP span a complete Hilbert space, and thus are conveniently described by Pauli 
matrices 

d\di + d\d 2 «-+ I = r° , (3.6) 
d\d x - d\d 2 t z , (3.7) 
d\d 2 + d^ 4-» t x , (3.8) 
d\d 2 - d\d x <-> ir y , (3.9) 

where I is the unit matrix. A suitable choice of a basis for the electron operator products c\ a Cj a — 1, 2) is: 

01 = c\ci - c\c 2 = c\c a + c\c e , 

2 = \ (c\ci - cJc 2 J = j (c\c - cjc e ) , 
Ei = c\ci + c\c 2 = c\c e + c\c a , 

E 2 = c\c 2 + c' 2 ci = c\c e - cj,c , (3.10) 

where Ei are symmetric and Oi antisymmetric under the transformation 1 2. 

Now the following assumptions are made: 
(i) the interaction between the atom and the electrons is invariant under the parity transformation between left and 
right. 

The energy of the atom at left and right positions may be different thus the tunneling atom may sit in an asym metr ic 
potential well. The assumption (i) is obvious for the screening inter actio n Caee the term proportional to V3 in Eq. ( |2.8j )) 
and the assisted tunneling exhibits the same invariance (see Eq. (2.21)). til This symmetry poses a constraint on the 
possible combinations of operators in the Hamiltonian: 

H z = t z (V 1 z 1 +V 2 z 2 ) , (3.11) 
H x = t x (V*E x + V 2 X E 2 ) , (3.12) 
H y = T y (V y 1 + V 2 V 2 ) , (3.13) 
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and the last combination, proportional to r° is not given since it can be incorporated into the potential scattering. 
The Vf's (i = 1, 2, a = x, y) are the appropriate couplings for assisted tunneling. The simplest form of the screening 
interaction is described by Vf, where the scattering on the atom is s-type. 

We note at this point that the most general form of the Hamiltonian may contain operators which can not be 



expressed simply by means of the operato rs in Eq. ( 3.1C ) because the operators c\ and C2 have simple specific 
structure in the momentum space (see Eqs. (2.4) and ( 2.22Q ), and interactions containing more elaborate momentum 



dependence can not be described in terms of these operators. 

A specific commutative model can be defined by adding to the general assumption (i) the followings 

(ii) The electron hole symmetry holds for the kinetic energy part of the Hamiltonian: g e {e) = g e (—e) = Q Q (e) = £? (— e). 

(iii) Additional symmetry: The interaction is also symmetric under the transformation: c\ *-* c|, C2 «-> c\. As discussed 
above this symmetry transformation is connected to the approximate local e-h symmetry of the Fermi surface. 

(iv) Time reversal symmetry, which connects the creation operators of the time reversal states A and A T , — > a\ T , 
and replaces any C-number by its complex conjugate as c — ► c*. 

Assuming that the wave functions of the tunneling atom are real (A T = A) 

dj = di , dl = (h ■ (3.14) 



Using the definition of time reversal symmetry for the free electron operators (ck) = c_k and Eq. (2.1) we find 
that (c±fc) T = c±k, and similarly c\ = c\ and = c%. 

(v) Moustakas and Fisher pointed out, that in the present case the transformation ci — > e'^ci and C2 — > e l< ^ 2 C2 
also leaves the Hamiltonian invariant. This symmetry is connected to the conservation of the number of electrons in 
channels 1 and 2. It actually follows from the previously mentioned assumptions (i — iv). 

For completeness we mention that in assumption (iii) the left-right transformation for the electron could be replaced 
by performing it in the atomic variable (d± <-> d%) instead of c\ «-* C2. In that case, however, the potential well must 
be symmetric, thus the atomic energies must be degenerate (A z = 0). For symmetry (iii) only the exchange of the 
relative positions of the HP and the electrons are important. 

As a consequence of the time reversal symmetry (iv) V% — 0, Kf = and V± = 0. The additional symmetry (iii) 
rules out all the terms except those proportional to 0±, thus that excludes the assisted tunneling. In that way the 
symmetries (iii) and (iv) allow only the interaction term of coupling Vf, thus they ensure the commutativity of the 
model even for an asymmetric TLS Hamiltonian (A 2 ^ 0). 

All the above consideration hold not only for the initial Hamiltonian, but for all the renormalized vertices. Thus the 
commutative model remains commutative during the scaling as a consequence of symmetries (ii)-(iv), and no terms 
can be generated which are different from and Ao- As also mentioned in the Introduction the elimination scheme 
of Moustakas and Fisher also violates this symmetry property of the Hamiltonian and generates terms even for the 
commutative TLS with the maximum symmetry (i)-(iv), which are not allowed. 



IV. INSTABILITY OF THE MARGINAL LINE DUE TO ELECTRON-HOLE SYMMETRY BREAKING 



It has been argued in Sec. II that the marginal line of fixed points is unstable due to breaking of the e-h symmetry. 
There is, however, an infrared cutoff appearing as the energy splitting Ao between the symmetric and antisymmetric 
heavy particle states as a consequence of the direct, spontaneous hopping between the left and right states. The 
asymmetry of the potential well A z further increases that cutoff to A = (Ag + A 2 Z )^ 2 . Thus the same direct hopping 
Ao is responsible for the instability of the fixed point and for blocking the scaling at lower energies. It will be shown, 
that for a TLS with realistic parameters in metals the scaling is restricted to the neighborhood of the marginal line 
which is far from the two channel fixed point. In order to demonstrate this statement in this section the notation 
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of the TLS literature is usedJS which results in much more symmetrical scaling equations. With these notations the 
interaction part of the Hamiltonian is 

H ** = E E^ c -</^t\ (4.1) 

i=0,x,y,z a,0 

where a, (3 = 1,2 label the orbital degrees of freedom and the spin index a has been restored. (Note that the matrix 
a 1 is acting in the orbital indices of the conduction electrons.) The V l, s denote the TLS-electron couplings, and a 
denotes the unit matrix. Note that now V° is initially set to zero, since its effect is incorporated in the local electronic 
propagators. Assuming g e (0) — Q o (0) = go for the sake of simplicity, the dimcnsionless values of the couplings are 
defined as v* = QqV 1 . This approximation is not essential, but makes the equations more transparent. These couplings 
of Ref. |^ can be expressed with the notations of the previous sections as V z — > V 3 , V x — > Ai and Vj., V% = 0. 

In the e-h symmetric case a e . Q = the scaling equations of the next to leading logarithmic order are&Eij 



dv° 

dx 
dv x 

dx 
dv y 

dx 
dv z 
dx 



= 0, 






= Av v v z 


-8v x ((v y ) 2 - 


- (v z ) 2 ) , 


= Av z v x 


-8v v {(v x ) 2 - 




= 4v x v v 


-8v z ((v x ) 2 ^ 


- K) 2 ) , 



(4.2) 
(4.3) 
(4.4) 
(4.5) 



where x = log Do /D with the scaled band width D, and the scaling of the splitting is described by 

^ = -8A [(^) 2 + (^) 2 ] . (4.6) 



(In Rcf. H instead of notation Ao the parameter A x = Aq is used.) Typical starting values of these parameters can 
be estimated as v z - 0.2, v x < 1CT 3 - 1CT 4 , v v = 0, A - 1 - 10 K, D - 10 5 K - 10 eV. 

The breaking of e-h symmetry can be taken into account in two different ways: (i) a slope in the density of states 
and (ii) different lower and upper cutoff, D up and D^ ovm . (One can also imagine a combination of the previous two 
cases). j— - 

First we consider the simplest case (i). In this case the general scaling equations are quite complicated,E2l however, 
in the limit \v z \ ^> \v x \, \v v \, they can be linearized in the 'small couplings', Aq, v x a nd v y , and one can convince 



himself easily that the discussed diagrams only modify the linearized version of Eq. (4.3). Then using the result given 



by Eq. (2.23) and combining that with Eqs. (4.2-4.4) one gets 

^A = -8K) 2 A , (4.7) 

^ = Av y v z - 8v x (v z ) 2 + 5a A v z2 , (4.8) 
dx 

rj v V 

— = Av x v z - 8v y {v z ) 2 , (4.9) 
dx 

where 6a — a e — a and v z rs const. These linear differential equations can be easily solved and one obtains for not 
too large v z 's: 
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Ao(») - A j , (4.10) 

z / n \ 4u*-8(d*) 2 

^(x)«fo^A ^J . (4.11) 

The scaling stops when (Aq/D)(x) ~ 1. Using the typical values given above and a ~ 1 (See App. |A|) the value of 
the generated coupling v x at the freezing out of the TLS can be estimated as 

< al ~ 10- 3 - 10- 4 , (4.12) 

which is too small to give an observable effect. This situation is sketched in Fig. ^. 

In the last part of this section we briefly discuss case (ii) with a ey0 = 0, f? e (0) = f? o (0) but with D up ^> Aow for 
both the even and odd channels, which means a very asymmetrical band. In this case a two-cutoff scaling can be 
used. The derivation of the scaling equations can easily generalized for the case of this asymmetrical cutoff. Assuming 
-D U p > Aow the integration now has to be carried out first for the energy range Aow < \u>\ < Ap and the couplings 
obtained by scaling must be then used for a symmetric model with D = A ow as starting values. Now the results 



of Sec. II can be generalized by keeping only the first part of integral ( 2.23 ). Thus the coupling v x at D = Aow 
(A P > Aow) is 

V*(D = Aow) ~ A (0)(T/*(0)) 2 (jj- - -L^ (MO) - Qo(0)) • (4.13) 

As for a realistic metallic case A P ~ Aow ~ 10 eV and £> e (0) ~ ^ o (0) ~ 1/D the corresponding dimensionless assisted 
tunneling can be estimated as v x (D = Aow) ~ 10 -6 . That initial value is too small to be increased essentially by 
scaling to obtain a value at the lower cutoff the scaling D = Aq which is comparable with the two channel fixed point 
v x w 0.25 . Indeed using Eq. (2.17b) of Ref. || one obtains 

v x {D = A ) = t-^Aow) cosh(4w z x) , (4.14) 

where x — log(Aow/Ao) ~ 10, thus 

v x (D = A )-6-10~ 4 , (4.15) 

which is still smaller by a factor ~ 10 -3 then the fixed point value v x ~ 0.25 . Finally it should be mentioned that the 
spontaneous tunneling rate Aq also renormalizes downward, but not more then 1-2 orders of magnitude, and since 



the typical renormalized value of Ao is ~ IK the estimation (4.15) is not essentially modified by the renormalization 
of A . 

In the discussion above it was assumed that Ap/Aow can not be very large. That is not the case-of a degenerate 
semiconductor, as a Pbi-^Ge^Te crystal with i<1, where orbital Kondo effect has been observed.EfJ In such a case 
even Ap/Aow ~ 100 can be realistic and also a can be present. In that case an essential strength of the assisted 
tunneling can be developed by scaling, even if it is negligible in the unrcnormalizcd Hamiltonian. That problem 
deserves further studies. 



V. CONCLUSIONS AND DISCUSSION OF THE PATH INTEGRAL TECHNIQUE 

In the previous sections the instability of the commutative TLS model against the breaking of the e-h symmetry 
has been studied using the multiplicative renormalization group. We found that without this symmetry breaking the 
line V x = V y = is marginally stable. 
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The instability of this marginal line has been shown by Moustakas and Fisher ,0 who pointed out that the original 
commutative TLS model has an artificial int ernal symmetry, connected to the conservation of the electron numbers 
in the left and right channels defined in Eq. ( [2.22 ). In the present paper we have shown that the model possesses an 
even more restrictive additional symmetry reflecting the role of e-h symmetry. This symmetry ensures the stability 
of the previously mentioned marginal line in the presence of e-h symmetry and in the absence of assisted tunneling. 
The breaking of the e-h symmetry combined with the screening V z and the spontaneous tunneling Ao leads to the 
generation of assisted tunneling which drives the system in the direction of the two-channel Kondo fixed point where 
the two channels are due to the double degeneracy of the conduction electrons because of their real spin. However, 
as it has been shown in Sec. IV, this fixed point can never be reached by the present mechanism as the driving force 
for the instability and the energy scale stopping the scaling are both provided by the spontaneous tunneling, Ao. 

It has also been pointed out that there is an ambiguity in the definition of the e-h symmetry transformation. Our 
symmetry (in) is valid for the kinetic part of the Hamiltonian in an arbitrary system nearby the Fermi surface if 
the curvature of the Fermi surface and the momentum dependence of the conduction electrons' energy are neglected, 
i.e., a local e-h symmetry is assumed at the Fermi surface, and the two minima of the TLS are close to each other: 
Rkf <C 1. As pointed out, this kind-of transformation can also be realized by an exact e-h symmetry transformation 
of a half-filled tight binding model.tll 

We have also estimated how far the system can get from the marginal line. For a conventional metal the dimen- 
sionless measure of the e-h symmetry breaking term can be estimated to be of the order of one, a ~ 1. Then the 
generated assisted tunneling at the freezing out of the TLS motion is of the order of v x ~ v y ~ 1CP 3 — 10~ 4 , which 
arc too small to produce an observable effect and the system is still not displaying the two-channel Kondo behavior. 
On the other hand, for degenerate semiconductors the high asymmetry in the lower and upper cutoffs_may generate 
a much larger v x via A , which might be responsible for the Kondo effect occurring in such materials r3 

The above discussion leads to the conclusion that for a TLS in an ordinary metal a strong enough initial assisted 
tunneling is required in addition to the Hamiltonian (2.5) to get a scaling into the intermediate vicinity of the two- 
channel Kondo fixed point and to have logarithmic minima in the electrical resistiyity.El The required assisted tunneling 
can also be generated via virtual excitations to the excited states of the HP.E3 The latter mechanism seems to be 
very promising, however, further studies are needed to check that these virtual hoppings do not generate a too large 
splitting Ao which might block the formation of the resistivity minimum. 

We stress at this point that the role of Ao is completely different in the commutative model and the non-commutative 
model with V x ,V y ^ 0.Q In the latter case A plays basically the role of an infrared cutoff parameter in the scaling 
and it drives the system from a non-fermi-liquid to a fermi-liquid transition^ This crossover has also been observed 
in recent point contact measurements, where it has also been possible to tune the crossover parameter ~ Ao by 
electron migration.EH This cross-over, of course, is only observable if Ao is not too large to stop the evolution of the 
non-Fermi liquid properties. We also mention that in the very strong coupling region Ao might become irrelevant, 
thereby playing again a qualitatively different role. This situation, however, is very unlikely to occur for realistic 
model parameters. nr - 1 

The present calculations support the general considerations of Moustakas and FisherE!rE3 concerning the instability 
of the marginal line, there are, however, essential differences. The source of these differences can be due to the 
ambiguities in the scaling procedures based on the path integral method. lSEJEj Since this technique is well established 
in the literature, we only briefhi-.sketch the procedure used. 

In the path integral methodEa the partition function of the system is written in an imaginary time path integral 
form, and is factorized as Z — Z\Z2, where Z\ involves only the screening operator V z , while all the hopping terms 
are pushed into the second term, Z2, which is calculated perturbatively in the -toppings. Following the same lines as 
Nozieres and de Dominicis in the solution of the X-ray absorbtion edge problemE3 one can integrate over the electronic 
degrees of freedom and obtain an effective functional for the TLS path, containing logarithmic interactions between 
the (spontaneous and assisted) flips of the TLS.C3 Finally, one can generate scaling equations for the 'fugacities' 
(corresponding to the different couplings and hopping amplitudes) by reducing the high energy cutoff intercoming in 
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the conduction electrons Green's functionSB 

The ambiguity of the path integral technique lies in the last step of the procedure, i.e. in the elimination of the 
high energy degrees of freedom, which corresponds to the change of a small time scale To ~ 1/ej? in the path integral 
scheme. 

In general two different elimination procedures can be applied: . .. . 

(i) One eliminates events of any kind which are closer in the imaginary time than the cutoff t rWA 

(ii) One eliminates only those events which are ckispr. to each other than to and in which electrons are involved.Ea 
Method (i) is followed by Moustakas suad FisherOO and in many other applications. It has been applied, e.g., 

for the problem of ID disordered metalstj where a special scheme had to be applied to avoid artificial generation 



of interactions. The aforementioned case of Ref. 21 is somewhat similar. On the other hand method (ii) rather 



corresponds to Anderson's poor map-, scaling, where the electronic heat bath is eliminated step by step via the 
reduction of the bandwidth cutoff D&3 

To show that method (i) delivers unphysical results let us consider the path of the TLS in Fig. |^ where a potential 
scattering occurs immediately after a spontaneous tunneling. Using method (i) Moustakas and Fisher replaced the 
event inside the box by a generated assisted tunneling interaction (see also Fig. 2 in Ref. |2l]). The artificial nature of 
this assisted tunneling is obvious since it is generated even if the dynamics of the TLS and the electrons are completely 
decoupled (v x = v y = v z = 0). This generated interaction, of course, can not be related to any measurable dynamical 
process. On the other hand, method (ii) does not generate this artificial interaction. 

The path integral method and the diagrammatic renormalization group are essentially different with respect to 
dynamical quantities. In the path integral approach, in principle, both methods (i) and (ii) can be followed, to 
calculate thermodynamical quantities even if method (i) is more frequently used. As a matter of fact in Appendix B we 
show that the scaling equations for the two-electron scattering obtained by methods (i) and (ii) are related by a formal 
transformation in the small coupling limit. The relationship between the different dynamical quantities is, however, 
much more complicated. To see this it is enough to consider the HP Green's function D(t) = (df (T)d 2 (T)dJ (0)<ii(0)) . 
While this Green function is invariant under the transformation (ii), using method (i) it acquires corrections of the 
type 6D ~ (df (T)d 2 (r)c+(T)c i (r)^(0)di(Q)) and - (df (T)d 2 (r)cT (0)cj(0)dj(0)di(0)) already in the first step of the 
renormalization group (i,j = 1,2), and under the subsequent steps it becomes a very complicated object composed 
from different many particle propagators. 

The multiplicative renormalization groop, on the other hand, guarantees the invariance of the different Green's 
functions. Thus the best way to choose the correct elimination procedure in the path integral technique is to compare 
it to the perturbative results. This comparison gives also the result that the it is method (ii) which gives the correct 
dynamical quantities. Therefore the scaling equations derived in Refs. [2l] and ^ must be interpreted with care. 

This conclusion may serve as a hint for other models where the strength of dynamical processes must be determined. 

Finally, we shortly discuss the case of two-electron scattering A 2 . It has been raised in Refs. ^l|,^2|that this process 
becomes relevant at the two-channel Kondo fixed point. While this statement seems to be correct in the highly 
anisotropical limit it is questionable in the strong coupling limit, where both the Bosonization approach (used in 
these works) and the Anderson- Yuval type scaling analysis loose their validity. On the other hand, this process has a 
small amplitude initially and is highly irrelevant in the weak coupling limit. This should be contrasted to the splitting 
of the TLS which is relevant in both the strong and weak coupling limits. Therefore, even if the two-electron scattering 
would be a relevant process at the Kondo fixed point, the dominant process driving away the TLS from the non- Fermi 
liquid fixed point will be the splitting for a realistic TLS, usually discussed in the literature. 

The generation of the A 2 interaction term by the e-h symmetry breaking can be analysed very similarly to that of 
the assisted tunneling process. Our investigations show that the generated A 2 is of the order of v z a 2 A and in the 
region v z < 0.4 the two-electron scattering is frozen out before it aquires an observable amplitude. On the other hand 
the role of A 2 in the strong screening case is still not completely clear and deserves further studies. 

To summarize, the general conclusion of the present study is, that the breaking of the e-h symmetry may have 
important consequences in cases, where the e-h symmetry ensures the cancellation of some diagrams (validity of 
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Ward identities) and thus the stability of some fixed points. In this way, the e-h symmetry breakiajprpay change the 
universality class of the problem, as it actually does i n ca se of the two-impurity Kondo problernElu3'E3E3 and in some 
other realizations of the two-channel Kondo problem. t2rllil 
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APPENDIX A: CALCULATION OF THE PARAMETERS a E ,o 



In the present Appendix we show that the e-h asymmetry parameters a e and a are generated even in the absence of 
potential scattering due to band structure effects. A simple free-electron like parabolic band with mass m is assumed 
for the electrons. The V3-term of the Hamiltonian can be derived by assuming a local electron-TLS interaction and 
it can be written as 



H 3 = V(g(R/2)- e(-R/2))r 2 



(Al) 



where g(r) = ip + (r)ip{r) denotes the density of the conduction electrons and we assumed that the average potential 
set by the TLS for the conduction electrons has already been taken into account in the calculation of the parabolic 
electronic dispersion. As the field operators are taken at the atomic sites a point-like s-wave TLS-clcctron interaction 
is assumed. In the conventional way an upper cutoff can be introduced by using a form factor in the momentum 
representation which also makes the interaction non-local. The interaction part (Al) can be rewritten by using the 
notations ip e ( ) = 3 (ip (R/2) ± ip (— R/2)) as 



where the field operators can be written in a momentum representation as 



(27T) 



(A2) 



(A3) 



Note that these operators are different from the c ej0 operators defined in Eq. (2.4) Then the spectral function for the 
Green's function G a b = — (T {ip^tpb}} = 8 a ,bG a , (a, b — e,o) can be expressed as 



(2tt) 3 



d 3 kS(LU 



2m' 



cos 2 (kR/2) 
sin 2 (kR/2) 



These integrals can be easily evaluated by taking the z-axis parallel to R, and one finally obtains 



Qe(o) M 



g{u) 



1± 



sin R(2mu) 
R(2rriLu) 



1/2 



1/2 



(A4) 



(A5) 
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where q{lj) = vr?l 2 2 1 ^ 2 (2n) 2 ui 1 ^ 2 is the density of states for the free ele ctron s for one spin direction, and the factor 
in the bracket is related to the normalization factor N e ^ a (k) given by Eq.( |2.3| ). In Ref. |l| these factors taken at the 
Fermi surface are incorporated in the definition of the coupling V3 as (g e (^F)Qo(^F)) 1 ^ 2 ~ (1 — (sinkpR / kpR) 2 ) 1 ^ 2 - 
There complete k-dependence can only be taken into account by a k-dependent coupling breaking the e-h symmetry. 



Finally, from expression (A5) the a parameters at the Fermi energy can easily be calculated as 

dlnge,o\ 1 lztcoskpR _i , 
duj J ep ~ 21±sink F R / k F R 6p 1 ' 

with k F = (2me^ ) 1 / 2 . The parameters a e>0 are of the order of e^ 1 and, e.g., for kpR <C 1 one obtains a e — a /3 — 
l/2e F \ 



APPENDIX B: CONNECTION BETWEEN THE TWO SCALING EQUATIONS OBTAINED BY THE 

TWO DIFFERENT PATH INTEGRAL METHODS 



In this Appendix we show that there exists a simple formal correspondence between the smalLcoupling limit of the 
scaling equations of Rcfs. |2l| and ^2] and those obtained by Vladar, Zimanyi and ZawadowskiEM This connection is 
based on the artificial mixing of the TLS and conduction electron degrees of freedom. 

Throughout this Appendix, for the sake of compactness, following Ref. || we write the Hamiltonian in the following 
form: 



fj,,cr,a,l3 



(Bl) 
(B2) 



where [i = 0, x, y, z and the indices a, [3 — 1,2 refer to the orbital degrees of freedom of the conduction electrons. In 
the special case considered in Sec. IV the couplings take the simple form V^g = ■cr^. For the conduction electrons 
a constant density of states go is used for both channels a — 1,2. 



Introducing the matrix notations go VTg 
tive renormalization group ardl 



u M the leading logarithmic scaling equations obtained by the multiplica- 



d ?Lo 
dx 

dx 



0, 



= 2ie l3k v,v 



.jUk , 



(B3) 



with e 1 ^ denotingthe Levi-Civita symbol. These equations coincide with the small coupling limit of the scaling 
equations of Ref. Eq obtained by using elimination scheme (ii) in Sec. V. 
The new couplings corresponding to Refs. El] and G2\ can be expressed as 



A* - 



-1/2 



A 1 



-1/2 



(B4) 



where a summation must be carried out over repeated indices. This transform, 
the high-energy region oj « D the unperturbed 'dimensionless' TLS propagatortu g 



on can be motivated as follows. In 
(°)(o;) = ojQ^{uj) can be written 



l(i 



as g(°^(D) w (1-AV/-D) -1 Transformation (B4) corresponds to mixing this factor (which is part of the TLS Green's 
function) into the v erte x in an artificial way. 

Writing out Eq. (B4) explicitly and plugging into Eq. (B3) we obtain the following scaling equations for the v^'s 



d|0 = 

dx ~ D ' 

dx - J - k ~ D 



(B5) 



Replacing the simple couplings of Sec. IV into these equations and putting A v — A z = one obtains exactly the 

We note at this point as well that terms proportional to ■ 



expanded version of the scaling equations in Ref. |22j. We note at this point as well that terms proportional to ~ A J / D 
give no contribution to the universal properties of the the model and can not be taken into account by a simple 
renormalization group procedure since they vanish in the scaling limit D — > oo, v^(D re {), A l (D re {) = est, D IC [ being 
a reference energy scale kept constant (it can be chosen to be the Kondo temperature, e.g., in the non-commutative 
model) . One can also show easily that the solution of the scaling equations derived from terms like ~ -jij In ^ does not 
repr odu ces the perturbative results and is inconsistent with the scaling hypothesis. Therefore the terms generated in 
Eq. (B5) by the transformation should vanish in the scaling limit. 



APPENDIX C: MIXING OF THE HIGH- END LOW-ENERGY STATES IN A SIMPLE POTENTIAL 

SCATTERING MODEL 

In this appendix we show how the original low- and high-energy degrees of freedom are mixed up by a simple 
potential scattering and how the elimination procedure during the renormalization group procedure should be carried 
out. For this purpose we investigate the simplest model possible: the scattering of spinless Fermions by a Dirac delta 
scatterer. In this case only s-wave scattering must be considered and the Hamiltonian can be written as 

H = H + Hu , 

#o = 5> lC +c„, (CI) 



a. in 



where we assume for the sake of simplicity that the spectrum of the unperturbed Hamiltonian H is given by e„ = 
^(n + 1/2), {n = —N,...,N}, L being the linear size of the system, and c+ denotes the creation operator of a 
conduction electron with label n. This assumption corresponds to a constant density of states (DOS) between the 
high- and low-energy cutoffs ±D = ±2irN/L. 



Eq. (CI) can be diagonalized easily via the unitary transformation 

c e = ^7(e,n) Cn , (C2) 
n 

where the operator c e annihilates an electron with energy e and 7(e, ri)"f(e', n) = 8 e>e i. The energy spectrum is 
determined by the self-consistent equation 



l = yE— ■ ( C3 ) 
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Then the factors 7 can be expressed as 



/ \ U 5(e) 

7(e.n) = y— L f-. (C4) 



where the constant S(e) is determined by the unitarity condition 



Eqs. (|C3|), (|C4|), and QC5J) constitute the complete solution of the problem. 

Our main purpose is to determine the exact eigenstates and the local density of states (LDOS) for the Hamiltonian 
above, which is given by the spectral function of the Green's function 

G(r) =-<Tc(0,t)c + (0,0)) , (C6) 

where c(0) = C^2 n c n )/^/L. This can be easily expressed in terms of the parameters j(e, n) as 

eH^%- £ )^l 7 ( e ,n)j . (C7) 

The last factor in this expression is nothing but the squared of the wave function amplitude A e — \ip e (0)\ 2 of the state 
with energy e. 

As one can convince himself easily from the investigation of Eq. ( CS ) the spectrum of the total Hamiltonian consists 
of 2N nearly equally spaced states forming a continuum with energies —D<t<D and a single bound stateo with 
energy 

on 

*B=±D± e±1/eQU _ i , (C8) 

where qq — l/2n and the '±' sign refers to attractive and repulsive interactions, respectively. This bound state gives 
a finite contribution to g(u>) with a weight 

'™ <0)|2 - w W) ' (C9) 



While the energy spectrum of the continuum is hardly affected by the potential scattering (see Eq. (|C3|)), the wave 
function of the states at the impurity site is, and therefore the continuous part of the LDOS becomes strongly energy- 
dependent. The selfconsistency equations for the states in the continuum can also be solved exactly by using the 
identity 

n— — 00 x / 



Then the self-consistency equation (C3) simplifies to 
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In the limit L, N — > oo, D = est the continuous part of the spectrum is given by 



2tt , . , 2 p 
e = -^(« + 1/2) + -<5„ , 

, c 2 1 £ + e 

ctg<3™ = 77 In , 

u tt D — e 

and the coefficients 7(e, n) can also be determined exactly: 

2 1 



7(e,n) 



£ /(e) e - £ r, 



with 



1 . D + e 2 
— In 

tt L> - e J7 



1 



1/2 



(C12) 
(C13) 

(C14) 
(C15) 



The wave function amplitude of the state with energy e can also be determined in this limit: 



A, = \M0)\ 2 = t 



U D + e V U 

— In 1 H 

2f D-6 / 4 



(C16) 



Replacing these amplitudes into Eq. ( |C7| ) we reproduce the results of the Green's function calculations Eq. ( [2.15 ) and 
one can easily calculate the parameter a as well: 



2U 



1 



tt(1 + U 2 /4) D 



(C17) 



The calculation above reveals the origin of the e-h symmetry breaking in the local properties. Eqs. ( |C4| ) and (|C16| ) 
tell us that the original high-energy states of Hq are mixed into the true low-energy excitations of the model with 
potential scattering. If one tries to do the renormalization group transformation by eliminate the high-energy states 
of Hq (D — > D') one immediately deforms the wave function of the true low-energy excitatio ns an d changes their 
coupling to the TLS drastically. This manifests in the change of their amplitude at the origin (C16) and the cutoff 
dependence of the parameter a (C17). Therefore the correct procedure is first to diagonalize the total Hamiltonian 
Hq + Hjj and then eliminate the high-energy states of the model. 



APPENDIX D: THE E-H SYMMETRY ON A CUBIC TIGHT BINDING MODEL 

The half-filled cubic tight binding band shows the e-h symmetry = — et+K with K = (7r/a, 7r/a, ir/a). In this 
cubic case the lattice can be divided into two sublattices A and B. The e-h transformation can be given in the 
site representation of the electrons' creation and annihilation operators as at — > af and bj — > —b~^ (icA, jeB). If 
one confines artificially the positions of the TLS onto the lattice sites even the non-commutative model is invariant 
under the above e-h transformation combined with a left-right exchange of the. |&Lectronic onsite operators similarly 
to Ref. This situation is similar to that of the two-impurity Kondo model. E30 

The situation is drastically changed if the tunneling atomic sites are not at the lattice points. For the sake of 
simplicity we assume that the the TLS sites ±R/2 are between two neighboring lattice points labeled by (±a/2, 0, 0) 
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and the corresponding operators are denoted by a + and b + , respectively. In this case the electronic field at ±R/2 is 
a combination of the atomic orbitals at sites corresponding to ±a/2. Assuming that the atomic orbitals ip are s-tyme 
the field operators at sites ±R/2 are 



*+(±R/r) = A 



\aTR\ 



aTR\ 



(Dl) 



A being a normalization constant. The interaction operators Eq. ( |2.8| ) and (2.21) can now be expressed in terms of 
v E'(±R/2) instead of cj.,2- If the overlap of the wave function -0 ^ l a T-R| ^ ^ / |-a=F-R| ^ j g no ^- ne gi ec ted then the assisted 



tunneling part of the Hamiltonian Eq. (2.21) is not invariant under the e-h transformation but it is invariant under the 



left-right exchange. Thus if the positions of the tunneling atom are not on the lattice sites then only the commutative 
mod el s hows the combined e-h invariance. This transformation is, however, defined differently from the one used in 
Sec. Ill and the difference is in the different momentum dependence of the transformation. We mention at this point 
that the cubic tight-binding model has a very large symmetry. Instead of the e-h transformation discussed in this 
Appendix one can also use an e-h transformation which connects electron, and hole states at the same part of the 
Fermi surface, the exact analog of the symmetry transformation in Fig. [|EjO 
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FIG. 1. The two first order vertex corrections for the heavy particle (dashed line) and conduction electron (solid line) 
interaction. In diagrams (i) and (ii) the intermediate state contains an electron and a hole, respectively. 



FIG. 2. Schematic representation of the interaction processes between the conduction electron and the HP. In the screening 
process (i) the atom sitting in one of the wells scatters the electron while in the assisted tunneling process (ii) the electron 
scattering induces the jump of the atom between the two wells. 



FIG. 3. Time ordered diagrams generating the assisted tunneling process in the presence of appropriate e-h symmetry 
breaking. The notations are identical to those of Fig. 1. The cross indicates a spontaneous tunneling between the two positions 
of the HP. The labels on the electron lines refer to the even and odd parity electron channels with respect to reflection through 
the center of the system. 



FIG. 4. (a) The electron and hole states (filled and open circles, respectively) connected by the approximate e-h transforma- 
tion through the Fermi surface element dS. (b) Schematic of the additional symmetry of the standard TLS model consisting 
of an e-h transformation and a reflection through the center of the TLS in the electronic degrees of freedom. 
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FIG. 5. Sketch of the scaling trajectories of the TLS. The appropriate e-h symmetry breaking drives the TLS away from 
the marginally stable fixed line V x = (Ai = 0) towards the two-channel Kondo fixed point. The scaling is stopped by 
the renormalized splitting, and the final ground state is a Fermi liquid. The freezing out of the TLS is indicated by a light 
continuous line. 



FIG. 6. The diagram resulting in the generation of artificial interaction terms in the elimination scheme (i). The heavy line 
represents the motion of the heavy fermion. The artificial interaction is generated by the elimination of a spontaneous tunneling 
event (Ao) and a potential scattering close to it (Vi.,2). The lines with arrows represent the conduction electrons. 



22 



FIGURES 




(M) 

FIG. 1. A. Zawadowski et al. /'Instability of the marginal commutative model of tunneling 
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